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A comparative density functional plus dynamical mean field theory study of the pseudocubic 
ruthenate materials CaRuOs and SrRuOs is presented. Phase diagrams are determined for both 
materials as a function of Hubbard repulsion U and Hund’s rule coupling J. Metallic and insulating 
phases are found, as are ferromagnetic and paramagnetic states. The locations of the relevant phase 
boundaries are determined. Based on the computed phase diagrams, Mott-dominated and Hund’s 
dominated regimes of strong correlation are distinguished. Comparison of calculated properties to 
experiments indicates that the actual materials are in the Hund’s coupling dominated region of the 
phase diagram so can be characterized as Hund’s metals, in common with other members of the 
ruthenate family. Comparison of the phase diagrams for the two materials reveals the role played 
by rotational and tilt (GdFeOs-type) distortions of the ideal perovskite structure. The presence of 
magnetism in SrRuOa and its absence in CaRuOs despite the larger mass and larger tilt/rotational 
distortion amplitude of CaRuOs can be understood in terms of density of states effects in the 
presence of strong Hund’s coupling. Comparison of the calculated low-T properties of CaRuOs to 
those of SrRuOs provides insight into the effects of magnetic order on the properties of a Hund’s 
metal. The study provides a simultaneous description of magnetism and correlations and explicates 
the roles played by band theory and Hubbard and Hund’s interactions. 
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I. INTRODUCTION 

The notion that the electronic properties of crystalline 
materials can to a large degree be understood in terms 
of the energy bands arising from the solution of the 
Schrodinger equation for a single electron in a periodic 
potential is fundamental to condensed matter physics and 
its applications. Electrons are charged and the inter¬ 
electron Coulomb interaction cannot be neglected. Den¬ 
sity functional theory (DFT), in essence a sophisticated 
mean field treatment of electron-electron interactions, 
provides a very good approximation to the interacting 
electron problem, enabling the theoretical description 
from first principles of many properties of many com¬ 
pounds. However DFT does not describe all electronic 
properties of all materials, and the cases where it fails 
can be taken to define the “strong correlation problem”. 

One dramatic example of strong correlations is the 
“Mott” insulator [1]: a material in which the correla¬ 
tions are so strong that they lead to insulating behavior 
in situations where band theory predicts a metal. Less 
extreme cases, where the interactions do not drive the 
material insulating but do lead to strong renormaliza¬ 
tion of electron velocity relative to band theory, to large 
and strongly temperature and frequency dependent elec¬ 
tron lifetimes, and to the occurrence of magnetic order, 
have been extensively documented [2]. 

Transition metal oxides (TMOs) [1, 2] play a particu¬ 


larly important role in the investigation of electronic cor¬ 
relations. In many TMO materials the transition metal d 
shells are partially filled. Interactions between electrons 
in the d orbitals of a transition metal ion are character¬ 
ized by a sizable effective Coulomb repulsion UeS that is 
close in magnitude to the bandwidth W of the d-derived 
bands and leads to the formation of spin and orbital de¬ 
grees of freedom. As a result the physics of TMO materi¬ 
als often differs sharply from the predictions of DFT and 
involves an intricate interplay of charge, spin and orbital 
degrees of freedom, which is furthermore highly sensitive 
to details of the crystal structure. 

These issues have been intensively studied in the con¬ 
text of TMOs which crystallize in variants of the ABO 3 
perovskite structure. In the ideal perovskite structure the 
B site ions lie on the vertices of a simple cubic lattice; 
each B site ion is octahedrally coordinated by oxygen. 
Few members of the family of ABO 3 transition metal 
oxides crystallize in the ideal cubic structure: in most 
materials a mismatch between the size of the A and B 
site ions (“tolerance factor” less than one) leads to a com¬ 
pressive strain on the BO 3 network. This strain is typ¬ 
ically accommodated by a rotational and tilt (GdFeOa- 
type) distortion of the BOq octahedra that diminishes 
the width of the d-derived bands and lowers the degen¬ 
eracy of the d multiplets. 

Particular attention has been given to materials in 
which the H-site ion is drawn from the first transition 
metal row of the periodic table so that the 3d shell of the 
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transition metal ion is partially occupied. In these ma¬ 
terials the key physics is the correlation-induced metal- 
insulator transition (often referred to as the ‘Mott’ tran¬ 
sition [1] although actual atomic-scale physics may be 
more involved [3]). The prevailing understanding [2, 4] 
is that in most of these materials the basic “correlation 
strength” is related to the proximity of the material to 
the Mott transition (but see Refs. 5-8 for the exceptional 
case of the nickelates) while the rotational and tilting 
distortions play a key role in determining this proximity. 
For example, SrVOs (nominal valence d^) crystallizes in 
the simple cubic structure and is a moderately correlated 
Fermi liquid [9]. In CaVOa, a small-amplitude GdFeOa 
distortion occurs; the material is still metallic but more 
correlated than SrVOa [9-11]. In LaTiOa the nominal d 
valence is also however a larger-amplitude GdFeOa 
distortion is present and the material is a Mott insu¬ 
lator. In isoelectronic YTiOa the distortion amplitude 
and the insulating gap are larger than in LaTiOs [2, 12]. 
The differences between Sr- and GaVOa or between La- 
and YTiOs may be attributed to different amplitudes 
of the GdFeOa distortion. Theoretical work [13] showed 
that the key physics is a lifting of the degeneracy of the 
transition metal t 2 g levels; this is important because the 
critical interaction strength needed to drive a Mott tran¬ 
sition depends strongly on orbital degeneracy, see, e.g. 
Ref. 14 and references therein. (The differences between 
the V-based and Ti-based materials arise in part from 
difference in GdFeOa distortion amplitude and in part 
from the difference in relative electronegativities of Ti 
and V [15, 16].) 

However, proximity to a Mott insulating state is not 
the only cause of correlated electron behavior. In heavy- 
fermion materials a lattice version of the Kondo effect can 
lead to enormous mass enhancements and other exotic 
physics [17]. In transition metal oxides with nominal 
valences from d^ to d® the Hund’s coupling can play a 
crucial role in producing very large renormalizations even 
for materials far from a Mott transition [18-23]. 

In this regard transition metal oxides where the tran¬ 
sition metal is drawn from the 4d series are of particular 
interest. Because 4d orbitals have a much greater spa¬ 
tial extent than 3d orbitals, the effective bandwidth is 
larger and the UeS is smaller, suggesting that the 4d ma¬ 
terials are in general likely to be farther from the Mott 
state than the 3d materials. Although many of the 4d 
series TMO are indeed itinerant metals, signatures of 
strong correlations, such as enhancement of the specific 
heat [24, 25], magnetic transitions [26] and evidence of 
other unusual electronic phases [27], are clearly present, 
in particular in the ruthenate family [23, 26]. Further, 
some members of the 4d series (for example Ca 2 Ru 04 ) 
have been identified as Mott insulators [28]. Thus in the 
4d-series transition metal oxides the issue of the relative 
importance of Mott and Hund’s correlations remains un¬ 
clear, as does the role of the GdFeOa distortions. 

Here, we explore these issues by focusing on the two 
of pseudocubic ruthenates: SrRuOa and CaRuOa . Both 


crystallize in GdFeOa-distorted versions of the ABO 3 
perovskite structure; with the distortion amplitude be¬ 
ing larger in CaRuOa than in SrRuOa . SrRuOa is fer- 
romagnetically ordered (a rather rare behaviour among 
4d TMOs) below a Gurie temperature Tc 160 K while 
CaRuOa remains paramagnetic to lowest temperatures. 
On the applied side, SrRuOa is a convenient electrode 
material, widely used as a substrate and magnetic ingre¬ 
dient in heterostructures and spin-valves [29-32]. Basic 
scientific questions remain open, including their degree of 
correlation, the origin of the apparently non-Fermi-liquid 
properties evident in the optical spectra [33, 34] and the 
reason for the magnetism, in particular why the appar¬ 
ently less strongly correlated material SrRuOa is mag¬ 
netic while the apparently more correlated CaRuOa is 
not. There is also fundamental interest in obtaining a 
better understanding of ruthenates in general, because 
insights gained in the study of the pseudocubic materials 
may shed light on the unconventional superconductivity 
of Sr 2 Ru 04 [26] and the metamagnetism and other phe¬ 
nomena observed in Sr 3 Ru 207 [27]. 

The question of the correct physical picture of the 
pseudocubic ruthenate perovskites (whether they should 
be regarded as weakly correlated itinerant metals or 
as strongly correlated systems) is the subject of con¬ 
troversy. On the experimental side, photoemission ex¬ 
periments [35] do not detect Hubbard sidebands, sug¬ 
gesting that the materials are not in proximity to a 
Mott transition. However, an earlier interpretation of 
the photoemission spectroscopy [36] indicated that siz¬ 
able renormalizations occur at low energies [37]. Optical 
spectroscopy [33, 34, 38, 39] indicates strong deviations 
from Fermi-liquid behavior, while transport experiments 
reveal very low Fermi liquid coherence scales (7K for 
SrRuOa [40] and 1.5 K for GaRuOa [39]) and large mass 
enhancements. 

On the theory side, early analyses [41, 42] of the elec¬ 
tronic structure based on spin-dependent density func¬ 
tional theory [local spin density approximation (LSDA) 
or spin-dependent generalized gradient approximation 
(GGA)] correctly describe many of the magnetic prop¬ 
erties. The ferromagnetism in SrRuOa was interpreted 
as the result of a Stoner instability, and the presence of 
magnetism in SrRuOa and its absence in CaRuOa was re¬ 
lated to the Fermi-level density of states, which is higher 
and more sharply peaked in SrRuOa than in GaRuOs . 
However the DFT calculations do not account for the 
low coherence scales and large mass renormalizations. 
A more recent comparative study of magnetism using a 
range of band theoretic techniques including the density 
functional plus U method concluded that U = 0 gives the 
best description of the experimentally observed transition 
temperatures [43]. Within LSDA, properties of SrRuOa 
and CaRuOa under strain were calculated [44, 45] and the 
predicted occurrence of ferromagnetism in CaRuOa un¬ 
der tensile strain was recently observed [46] . On the other 
hand, many theoretical papers including the LSDA-l-17 
work of Rondinelli and collaborators [47] and several dy- 
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namical mean-field theory (DMFT) investigations [48- 
51] assert that correlations beyond LSDA/GGA are im¬ 
portant. 

The existing literature thus suggests that the challenge 
presented by the perovskite ruthenates is to develop a 
theory that includes the electronic correlations that pro¬ 
vide the experimentally indicated mass enhancement and 
other renormalizations without spoiling the good account 
of the magnetic phase diagram found in density func¬ 
tional calculations. In this paper we address this chal¬ 
lenge by performing a systematic density functional plus 
dynamical mean field theory study that includes realis¬ 
tic electronic structure and investigates a wide range of 
potentially relevant interaction parameters. We calcu¬ 
late the phase diagram in the (U, J) plane and by com¬ 
paring calculated and measured properties we locate the 
perovskite ruthenates in the correlated Hund’s metal re¬ 
gion of the phase diagram. Effective masses are found 
to be large and coherence scales small in the paramag¬ 
netic phase. The greater tendency to magnetic ordering 
in SrRuOs than in CaRuOa is accounted for in a manner 
similar to that found in earlier electronic structure calcu¬ 
lations. Our results also provide insight into the general 
issue of the effect of GdFeOa distortions on the effective 
correlation strength of Hund’s metal. The low frequency 
properties are affected by the density of states (which in 
the ruthenates is reduced by GdFeOs distortions) while 
the global and higher frequency correlation strength is 
controlled by the inverse bandwidth (which is increased 
by GdFeOa distortions). 

The rest of this paper is organized as follows. Section II 
describes the methods we use. Section III presents our 
computed ferromagnetic/paramagnetic (FM/PM) and 
metal-insulator phase diagrams and discusses the physics 
behind them. Section IV discusses the differences of two 
ruthenates in the mass enhancement and the self-energy 
and uses this information to locate the materials on the 
phase diagram of Sec. III. The magnetic phase of SrRuOa 
is analyzed in detail in Sec. V. Section VI presents a sum¬ 
mary and prospects for future work. Appendices provide 
details of the calculations 


II. CRYSTAL STRUCTURE, ELECTRONIC 
STRUCTURE AND MODEL 

CaRuOa and SrRuOa crystallize in a Pnma symme¬ 
try crystal structure related to the ideal cubic perovskite 
structure by a GdFeOs distortion corresponding to a tilt 
and rotation of each RuOg octahedron. The tilts and 
rotations alternate in a four-sublattice pattern. SrRuOa 
has a Ru-O-Ru bond angle of about 163°, in between 
the ideal perovskite Ru-O-Ru bond angle of 180° and 
the Ru-O-Ru bond angle of 150° observed in GaRuOa 
[52, 53]. 

Valence counting implies that in CaRuOa and SrRuOa 
the Ru is in the electron configuration with four elec¬ 
trons in the Ru 4d shell. The octahedral ligand held 


pushes the Cg levels up in energy so the relevant near- 
Fermi-surface bands are derived from Ru t2g symmetry d 
states with some admixture of the O 2p states. Because 
Ru is a second-row transition metal ion, the d states are 
expected to be more extended and the on-site interaction 
U is expected to be weaker than for the hrst-row transi¬ 
tion metal ions, indicating [23] that the materials are not 
in the charge transfer regime of the Zaanen-Sawatzky- 
Allen phase diagram [3]. We therefore adopt the “fron¬ 
tier orbital” approach in which the low energy electronic 
properties are obtained from a multiband Hubbard model 
with hybridizations and level splittings obtained from the 
near-Fermi-surface transition metal d-derived bands. 

The Hamiltonian takes the general form 

dd — dIkin 4" ddonsite? (1) 

where idkin describes the dispersion of the bands derived 
from the frontier orbitals and idonsite the additional in¬ 
teractions. The chemical potential is set to ensure that 
these bands contain four electrons per Ru. 

To define the near-Fermi-surface bands of idkin pre¬ 
cisely we use the non-spin-polarized generalized gradient 
approximation (GGA) as implemented in the Quantum 
ESPRESSO density functional code [54, 55] to obtain 
electronic band structures and then project the result¬ 
ing bands onto maximally-localized Wannier functions 
(MLWF) [56, 57] using the wannier90 code [58] (details 
are given in Appendix A). 

In most of this paper we construct idkin by projecting 
the Kohn-Sham Hamiltonian onto t 2 g-symmetry Wan¬ 
nier functions centered on the Ru sites. This procedure 
captures correctly all of the electronically active frontier 
orbitals and provides a reliable description of the phase 
diagram and quasiparticle properties. However, as will 
be discussed in detail in Sec. V, this procedure leads to 
an overestimate of the magnetic moment in the magnet¬ 
ically ordered phase. Obtaining a correct estimate of the 
ordered moment requires inclusion of bands derived from 
Ru Cg states. In our analysis of the magnetic state the 
Bg-derived bands are therefore retained, but because the 
Bg-derived states are far from the Eermi level, they are 
treated by the mean-field approximation used in Ref. 59. 

The t 2 g orbitals are treated dynamically. As usual in 
studies of transition metal oxides, the interaction Hamil¬ 
tonian is taken to be site-local and to have the rotation- 
ally invariant Slater-Kanamori form [2] . We use the form 
appropriate [23] for intra-t 2 g orbitals, since these are the 
primary focus of this work 

dionsite — U E 

OL 0(^0 

+ {U- 3J) E 
4- E 4- 

a//3 

where a, /3 are orbital indices and a is the spin index. Dif¬ 
ferent values of U and J have been used for ruthenates 
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in the literature. For Sr 2 Ru 04 , constrained LDA [60] 
gives {U,J) values of (3.1 eV, 0.7 eV) while the values 
(2.3 eV, 0.25 eV) [20] and (2.6 eV, 0.26 eV) [61] have been 
obtained from constrained random phase approximation 
(cRPA) method. As discussed above, for the perovskite 
ruthenates U values ranging from zero to rather large 
numbers have been employed. For this reason and be¬ 
cause the behavior of the model for general parameters 
is of theoretical interest, we consider a range of values for 
U and J in this paper. However we restrict attention to 
the regime U > 3J where the effective on-site interaction 
is positive in all channels. 

To treat the onsite interaction [Eq. (2)], we em¬ 
ploy single-site dynamical mean-field theory (DMFT) 
[62]. This method allows us to map the Hamiltonian 
[Eq. (1)] into a multiorbital impurity model embedded 
in a fermion bath. The impurity model is solved by us¬ 
ing the hybridization expansion version of the continuous 
time quantum Monte Carlo (CT-HYB) [63] implemented 
in the TRIQS code [64] for rotationally invariant inter¬ 
action using conserved quantities [65] to speed up the 
calculations. 

Care is required in the definition of the impurity model. 
Each Ru ion is at the center of an octahedron defined by 
six oxygen ions. The Pnma structure means that the 
local symmetry axes of a given RuOg octahedron are not 
parallel to the axes that define the global crystal struc¬ 
ture. If the Cg and t 2 g combinations of the d-derived 
states in a single octahedron are defined with respect to 
the global, rather than the local, symmetry axes, the im¬ 
purity model will contain off-diagonal terms which mix 
the different orbitals at the single-particle level. This 
causes a severe sign problem for the CT-HYB solver 
[66]. It is preferable to avoid this complication, follow¬ 
ing Ref. 67, by using a local basis with symmetry axes 
aligned along the octahedral directions appropriate to a 
given Ru, instead of the global axes. When restricted to 
the t 2 g manifold only, the MLWF approach used here pro¬ 
duces orbitals that are already aligned with respect to the 
local axes of the appropriate octahedron so the hybridiza¬ 
tion function is essentially diagonal. Thus if only t 2 g 
orbitals are retained the DMFT calculation is straight¬ 
forward: the impurity model with hybridization function 
defined directly from the projection of the Kohn-Sham 
Hamiltonian onto the Wannier states is solved for one 
Ru site. The self-energies for the other Ru sites are then 
constructed by applying appropriate rotation operators. 
If the MLWF procedure is applied to the full d manifold 
(including both t 2 g and Cg orbitals) then the resulting 
orbitals turn out not to be aligned to the local symmetry 
axes and an additional change of basis is required before 
solving the impurity model (see Appendix A). From the 
solution of the dynamical mean field equations we de¬ 
termine the phase (metal versus insulator, paramagnetic 
versus ferromagnetic) and some properties of the phases, 
in particular the quasiparticle mass enhancement and the 
magnetic moment. Details of our procedure for determin¬ 
ing different phases are given in Appendix B. 


III. QUALITATIVE PHYSICS 
A. Electronic structure: Density of states 



FIG. 1. Orbitally-resolved density of states generated from 
MLWF fittings of the DFT calculations for CaRuOs (SrRuOs) 
on the positive (negative) half plane. Each panel is corre¬ 
sponding to each of the t 2 g orbitals. The (xy, yz, zx) occupan¬ 
cies are (1.40,1.36,1.23) for CaRuOa , and (1.34,1.26,1.40) 
for SrRuOa , respectively. In the same order, the DOS at 
the Fermi level vf (states/eV per Ru atom) for each of 
the three t 2 g orbitals are (0.83,0.73,0.79) for CaRuOa , and 
(1.00,0.77,1.19) for SrRuOa . The vertical dashed line marks 
the Fermi level. 

Figure 1 shows the orbitally-resolved near-Fermi- 
surface density of states (DOS). The DOS of the two 
materials are similar, as expected from the essentially 
identical quantum chemistry, but the difference in the 
magnitude of the GdFeOa distortion occurring in the two 
materials leads to two important differences in the DOS. 
First, the t 2 g-derived bands in SrRuOa are approximately 
10% wider than those of CaRuOa (SrRuOa bandwidth 
« 3.0eV, as compared to « 2.7eV for CaRuOa). To 
the extent that correlation effects scale as the ratio of 
an interaction strength to a bandwidth, this would sug¬ 
gest that CaRuOa would be the more strongly correlated 
material. 

However, SrRuOa has the larger density of states at 
the Fermi level. This can be traced back to the van Hove 
singularity of the undistorted cubic structure, which hap¬ 
pens to lie very close to the Fermi level. Because the 
GdFeOa distortion lifts the degeneracy of the t 2 g lev¬ 
els, it splits the van Hove peak into three features. In 
SrRuOa the splitting is small and the density of states 
remains large. In CaRuOa the splitting is larger, leading 
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to a smaller Fermi-level DOS. To the extent that correla¬ 
tion effects are related to the Fermi-level density of states, 
this suggests that SrRuOa would be the more strongly 
correlated material. In particular, the Stoner model of 
ferromagnetism relates the presence of magnetic order 
to the value of the product of an interaction constant 
and the Fermi-level density of states [68], so the density 
of states difference would suggest (in agreement with ex¬ 
periment and with the DFT work of Refs. 41, 42, and 69) 
that SrRuOa is more likely to be magnetic than CaRuOa . 
Further, particularly in CaRuOa, the splitting creates a 
density of states peak below the Fermi level. The consid¬ 
erations of Ref. 67 building on previous work of Vollhardt 
and collaborators [70-74] suggest that this peak is unfa¬ 
vorable to magnetism. 


B. Phase diagrams 

Figure 2 displays the phase diagrams in the plane of the 
Hubbard U and the Hund’s coupling J, determined by 
the procedure described in Appendix B. We considered 
metallic and insulating phases, paramagnetic and ferro¬ 
magnetic phases. Antiferromagnetism was not studied. 

Focus first on the upper right panel, which presents re¬ 
sults for CaRuOs. We see that as the interaction strength 
U is increased at fixed J, there is a phase transition from 
a metal (which may be paramagnetic or ferromagnetic) 
to a Mott insulator. As J is increased at fixed t/, a tran¬ 
sition to a ferromagnetic metal phase occurs but no fer¬ 
romagnetism was found in insulating phases in the U and 
J ranges that we studied. From this phase diagram one 
can identify two regimes: at large values of U, near the 
metal-insulator phase boundary, properties are most sen¬ 
sitive to the value of the effective Hubbard interaction. 
Note in particular that at large J the phase boundary 
becomes a straight line with slope U — 3 J. The quantity 
U — 3J is the effective Hubbard interaction (correlation 
strength) relevant to the Mott transition because it gives 
the lowest energy cost for a valence change from 2d^ to 
high spin (Pd^. In contrast, far from the metal-insulator 
phase boundary, the Hund’s coupling J is the key param¬ 
eter: by increasing J a ferromagnetic phase is induced 
and there is a range of U in which the location of the 
critical boundary is only weakly dependent on U. 

The upper left panel displays our results for SrRuOa. 
The same phases are found but the difference in GdFeOa 
distortion amplitude causes the location of the phase 
boundaries to be different. To highlight the differences 
between the two calculations we present in Fig. 2(c) a su¬ 
perposition of the two phase diagrams. At larger U the 
phase boundaries are parallel (and controlled hy U — 3 J), 
with CaRuOa requiring a slightly smaller value of U to be 
driven into the Mott phase, as expected from the smaller 
bandwidth and larger t 2 g-level splitting following from 
the larger distortion in CaRuOa . However, in this region 
of the phase diagram the magnetic phase boundaries are 
very similar, and an extreme fine-tuning of {U, J) would 



FIG. 2. (Color online) Ferromagnetic metal/paramagnetic 
metal (FM/PM) and metal/Mott insulator (MIT) phase di¬ 
agrams. The dotted lines indicate the metal-insulator phase 
boundary, with the region above the lines being insulating. 
The heavy solid line (red on-line) indicates the boundary of 
the ferromagnetic region, while the light line (black on-line) 
separates the physically relevant positive effective interaction 
{U > 3J) from the unphysical negative interaction {U < 3J) 
region. (a,b) Phase diagrams for SrRuOs (a) and CaRuOs (b) 
obtained from the DFT-I-DMFT procedure described in the 
text. Open circles (black on-line) indicate ([/, J) points for 
which properties were computed and a paramagnetic metallic 
state was found. Closed circles (black on-line) indicate para¬ 
magnetic insulating solutions. Open squares (red on-line) in¬ 
dicate (U, J) points for which a ferromagnetic metallic state 
was found, (c) The phase diagrams for SrRuOa and CaRuOs 
plotted together. Ferromagnetic phase boundary of CaRuOs 
indicated by heavy dashed line (blue on-line), (d) Phase di¬ 
agrams for SrRuOs and CaRuOs plotted together, but with 
the phase boundary of SrRuOs rescaled by the ratio of the 
SrRuOs bandwidth to the CaRuOs bandwidth. Note that 
there is no FM insulating phase found in these phase dia¬ 
grams. 

be required to account for the fact that CaRuOs is a para¬ 
magnetic metal, while SrRuOs is a ferromagnetic metal. 

A much more significant difference between phase 
boundaries for the two structures is found at smaller val¬ 
ues of C/ < 3 eV. There, the CaRuOs phase diagram has 
a significantly smaller region of ferromagnetism than the 
SrRuOs one. Thus, in this regime much less fine-tuning 
of the parameters is needed to correctly account for the 
difference in magnetic behavior of the two materials. As 
discussed in Sec. IV, there are also other experimental 
indications suggesting that these two materials should 
be viewed as being in this Hund’s coupling dominated 
regime. We also note that the cRPA values of U, J found 
for the related Sr 2 Ru 04 material [61] are in the mod¬ 
erate I/, larger J region where CaRuOs is paramagnetic 
but SrRuOs is ferromagnetic. 

We now discuss further the qualitative origin of the 
observed differences between the phase diagrams of the 
two structures, by presenting in Fig. 2(d) a superposition 
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of the phase diagrams of the two materials, but with the 
U and J values for SrRuOa rescaled by the ratio r = 1.11 
of the SrRuOa to the CaRuOs bandwidth. At larger U 
the phase diagrams for both the metal-insulator transi¬ 
tion and magnetism coincide in the rescaled plot. This 
indicates that in the Mott-dominated region the physics 
is controlled by the ratio of the interaction strengths to 
the bandwidths and depends only weakly on for example 
the Fermi-level DOS. However, in the smaller U regime, 
the magnetic phase diagrams do not coincide, indicating 
that in this regime the physics is clearly not controlled 
solely by the difference in bandwidths. 

Instead, the substantial difference in the critical J re¬ 
quired to drive the ferromagnetic transition is associated 
to the DOS in the near-Fermi-level region. One impor¬ 
tant property is the value of the DOS at the Fermi level 
vp. In the standard Stoner theory [68] magnetism is asso¬ 
ciated with a value greater than unity of a dimensionless 
interaction parameter I constructed as the product of an 
appropriate interaction energy and the Fermi level den¬ 
sity of states (/ = Uvp > 1). Clearly the larger DOS 
in SrRuOa makes it easier for the Stoner parameter to 
exceed the critical value and as discussed by Mazin and 
collaborators [41, 42] spin dependent DFT calculations 
indeed indicate a Stoner parameter slightly greater than 
unity for SrRuOa and slightly less for CaRuOa . It is 
also worth noting that the Stoner theory is in essence 
a Hartree approximation. When the correlation is fully 
treated, other factors such as the energy derivative of 
the DOS at the Fermi level and indeed the structure of 
the DOS far from the Fermi level are also important and 
provide significant corrections to the simple Stoner esti¬ 
mate. Dynamical mean field studies of related models 
[67, 70-72] indicate that for systems with carrier concen¬ 
tration such that the d shells are less than half occupied 
such as Lai-xSrajVOs {d valence ferromagnetism 

is favored if the DOS peak is at or below the Fermi level 
[67]. A particle-hole transformation allows us to relate 
the results of Ref. 67 (derived for a system with valence 
near d^) to the ruthenates (valence d'^), concluding that 
in the ruthenates a DOS peak at or above the Fermi level 
favors ferromagnetism. Therefore, as seen in Fig. 1, the 
larger distortion of CaRuOa produces below-Fermi-level 
density of states peaks, thus disfavoring ferromagnetism, 
while in SrRuOa the DOS peaks concentrate at the Fermi 
level and ferromagnetism is favored. 


IV. SELF-ENERGY, MASS ENHANCEMENT 
AND QUASIPARTICLE LIFETIME 

In this section we study the electron self-energy and 
quasiparticle properties, choosing interaction parameters 
U = 2.3 and 3 eV as representatives of the Hund’s metal 
and Mott-dominated regimes respectively. These values 
are similar to those obtained for the related material 
Sr 2 Ru 04 from ab initio estimations using constrained 
DFT [60] and constrained RPA [20] approaches, respec- 



“n(eV) 



FIG. 3. (Color online) The imaginary part of the Matsub- 
ara self-energy ImE(ia;„) averaged over the three t 2 g orbitals, 
calculated at [/ = 2.3 (a) and 3eV (b) with J = 0.4eV, 
T — 0.0025 eV for SrRuOa and CaRuOs . The inset provides 
an expanded version of the low frequency regime, allowing 
the crossing point u)c to be distinguished. Calculations are 
restricted to the paramagnetic order. 


tively. We fix J = 0.4 eV as a representative value for 
which the Ca material is paramagnetic and the Sr mate¬ 
rial is ferromagnetic. 

Figure 3 displays the imaginary parts of the com¬ 
puted Matsubara-axis DMFT self-energies restricted to 
the paramagnetic phase over a wide range of Matsub- 
ara frequencies. We see that at low frequencies the self¬ 
energy of the Sr-compound is larger, indicating that for 
quantities dominated by low energies the Sr material may 
be viewed as more correlated. On the other hand, above a 
frequency ujc the self-energy of the Ca material is larger, 
reflecting the effect of the difference of bandwidths on 
the effective correlation strength at higher energies. The 
Hubbard U and Hund’s coupling J compete in this re¬ 
spect. At smaller U [Fig. 3(a)], the Hund’s coupling effect 
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FIG. 4. (Color online) The mass enhancement Z~^ com¬ 
puted as described in the text for SrRuOa (squares, black 
on-line) and CaRuOa (triangles, red on-line) at J = 0.4 eV 
and U = 2.3 eV (solid lines) and 3eV (dashed lines) as a 
function of temperature. In these calculations magnetism 
is suppressed: all results are for the paramagnetic phase of 
the model. Inset: comparison of mass enhancement (squares, 
green on-line, left y axis) and zero Matsubara frequency spin- 
spin correlation function x(j^^n=o) = Jq {Sz{r)Sz{0))dT with 
Sz = — NI°^) (circles, blue on-line, right y axis) com¬ 

puted for CaRuOa at temperature T = 2.5 meV ~ 29K with 
U = 2.3 eV and plotted as a function of J. 

is stronger, resulting in a wider range of low frequency 
(larger ujc ~ 0.4 eV) in which SrRuOa is more correlated. 
For the larger U close to the MIT phase boundary, Mott 
physics associated with U becomes stronger, as signaled 
by the decrease of ojc by a factor of ^ 4 [wc ~ 0.1 eV in 
Figure 3(b)]. 

We now turn to two key physical quantities charac¬ 
terizing quasiparticles in the metallic state, namely the 
effective mass enhancement m*/m (directly related in 
the local DMFT approach to the quasiparticle weight Z) 
and the quasiparticle scattering rate (inverse lifetime) F. 
These are defined from the real and imaginary part of 
the retarded self-energy on the real-frequency axis by: 

— = ^ = l-^ReS(cc + *0+)U=o, (3) 

m Z OLJ 

Inferring real-axis quantities from Matsubara-frequency 
data in general requires analytical continuation. If, how¬ 
ever, the low frequency properties are reasonably well 
described by the Fermi liquid fixed point (as is the case 
for the parameters we study), the low frequency limit 
of the real frequency self-energy may be inferred with 
reasonable accuracy from the data at small Matsubara 
frequencies, with I — Z~^ k. dhnS{iuJn)/dLOn\uj„^o 
Z~^r = —ImT, (iujn —^0). In practice, we extract Z 
and F by fitting a fourth-order polynomial to the first six 
Matsubara-axis data points for ImE(ia;„) and computing 
the needed quantities from the fitting function. 


Figure 4 shows the estimated mass enhancement for 
the two materials at the two values of U under considera¬ 
tion. The calculations are restricted to the paramagnetic 
state and show a strong temperature dependence, which 
is a manifestation of the low quasiparticle coherence scale 
associated with the formation and slow fluctuation of a 
local moment with 5" ^ 1/2 (in other words, with Hund’s 
metal physics) [19, 23]. As will be seen below, in the fer¬ 
romagnetic state of SrRuOa the temperature dependence 
is cut off because the ferromagnetic order quenches the 
slow spin fluctuations. To reinforce this point we show in 
the inset that the mass enhancement and zero Matsubara 
frequency impurity model spin correlation function have 
identical J dependence. Comparison to Fig. 1 shows that 
the local spin fluctuations involve an energy scale that is 
much lower than that of the bare density of states. 

In our simulations we were unable to reach tempera¬ 
tures low enough to observe the saturation of the mass to 
its T —>■ 0 limit in the paramagnetic phase. Nevertheless, 
we find the quasiparticles are becoming well defined at 
the lowest temperatures reached in our simulations. Fig¬ 
ure 5 presents the temperature dependence of F/T calcu¬ 
lated for CaRuOa at 17 = 2.3 and 3eV with J = 0.4 eV. 
Below T = 70K for [/ = 2.3 eV and below T = 30K 
for U = 3eV the scattering rate becomes smaller than 
temperature, which is indicative of coherence. We ex¬ 
pect that as the temperature is lowered further below the 
coherence scale the mass will saturate. From these con¬ 
siderations we estimate the T — 0 mass enhancements to 
be about 7 and 11 for the U = 2.3 and 3eV, respectively. 



T (K) T (K) 

FIG. 5. (Golor online) Inverse of the quasiparticle lifetime, F 
(Eq. (3)) calculated for GaRuOa at J = 0.4eV and 17-values 
indicated, divided by temperature and plotted as a function 
of temperature, for each of the t 2 g orbitals. Linearity cor¬ 
responds to Fermi liquid behaviour F oc The horizontal 
dashed line marks the boundary (F < T) of the coherence re¬ 
gion. The three orbitals differ slightly in occupancy; the closer 
the orbital is to half filling the more strongly it is correlated. 
[20]. The convention for the orbitals is given in Appendix A. 

We remark that Fig. 5 implies that the quasiparticle 
scattering rate F varies as up to T « 150 K, even 
though the characteristic Fermi liquid signatures in phys¬ 
ical observables (temperature independent mass/specific 
heat coefficient and magnetic susceptibility along with 



























J (eV) 

FIG. 6. (Color online) The mass enhancement of (a): SrRnOa 
and (b): CaRuOa vs. the Hund’s conpling J calculated at 
U = 2.3 eV for several different temperatures. For SrRuOs, 
the larger square indicates the J value closest to the onset of 
ferromagnetism at the corresponding temperature. 


quadratic transport scattering rate) are only evident be¬ 
low much lower temperatures (lower than the lowest tem¬ 
perature accessible in our U = 2.3, J = 0.4 eV calcu¬ 
lations). This behavior qualihes CaRuOa as a ‘hidden 
Fermi liquid’ [75, 76] in which although the temperature 
dependence of e.g. the resistivity deviates from above 
a very low temperature, the quasiparticle scattering rate 
remains ~ up to much higher temperatures and the 
deviation from the Fermi liquid temperature dependence 
expected for the resistivity is attributable to a tempera¬ 
ture dependence of the quasiparticle weight. 

We now use the experimentally measured low temper¬ 
ature specihc heat coefficient to help constrain the pa¬ 
rameter values. The measured low-T specific heat of 
CaRuOais « 74mJ/molK^ [77], implying a mass en¬ 
hancement m*I'm « 7 with respect to the DFT value. 
We note that the experimental value has a contribution 
(of unknown magnitude) from the electron-phonon inter¬ 
action, so should be regarded as an upper bound on the 
electronic contribution to the mass enhancement. Fig¬ 
ure 4 shows that has a marked dependence on U 



FIG. 7. Predicted orbitally-resolved spectral function of 
GaRuOs. We used U = 2.3 eV, J = 0.35 eV, and temperatnre 
T = 2.5 meV « 29 K. Dashed line is the orbitally-averaged 
spectral function for GaVOs with data extracted from Ref. 78. 


while Fig. 6 shows that Z~^ depends even more strongly 
on Hund’s coupling J. As the interaction parameters 
are not likely to change significantly between the two 
compounds, we assume that SrRuOa and CaRuOa are de¬ 
scribed by the same (t7, J) values. Requiring that the 
calculated mass enhancement for CaRuOa is close to but 
not higher than the measured mass, and at the same 
time that J be such that SrRuOa is ferromagnetic and 
CaRuOa is paramagnetic allows us to locate the materi¬ 
als on the phase diagram. 

At C/ = 3eV the phase diagram of Fig. 2 requires 
that 0.3 eV ^ J ^ 0.4 eV while for J in this range the 
masses resulting from our calculations are clearly above 
the experimental value (« 10 at T = 50K and clearly 
increasing as T is decreased; see also the lower panel of 
Fig. 8). Thus we argue that the combination of the mass 
and phase diagram are inconsistent with the possibility 
that the perovskite ruthenates are in the Mott-dominated 
regime. 

On the other hand, our results at smaller U indicate 
that the Hund’s metal regime can provide a good de¬ 
scription of the basic physics. A relatively wide range of 
J is found for which SrRuOa is magnetic and CaRuOa is 
not, while the sensitive dependence of the mass enhance¬ 
ment on J, with masses ranging from much smaller than, 
to rather larger than, the measured value, means that 
a reasonable J ^ 0.3 eV and U ^ 2 — 2.5eV (close to 
those found from constrained DFT and cRPA methods 
[20, 60, 61]) will account for the basic physics. 

It is interesting to look also at the correlated DOS. 
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We obtained the spectra at real frequencies using the 
maximum entropy method [79, 80] to continue the self¬ 
energies [81] and then using the continued self-energies to 
construct the orbitally-resolved spectra. The results are 
presented in Fig. 7 (solid lines). The densities of states 
of the three orbitals are quite similar. Comparison of the 
near-Fermi-level fine structure in the DOS to that shown 
for the noninteracting model in Fig. 1 reveals that the 
peaks seen at ~ —0.4 eV and ~ -1-0.35 eV in the xy and 
yz bands of the CaRuOa DFT DOS are renormalized dif¬ 
ferently. The negative frequency peak is pulled up much 
closer to the Fermi level, appearing at about —0.05 eV 
in the correlated band structure. The positive frequency 
peak in the correlated band structure occurs at approx¬ 
imately the -1-0.2 eV energy of the DFT band structure, 
thus is renormalized by a much smaller amount. The dif¬ 
ference is a manifestation of a large particle-hole asym¬ 
metry in the self-energy, whose influence on the spec¬ 
troscopy was discussed in [82] . The effect is a fingerprint 
of Hund’s metal physics, as recently noted in Ref. 83. See 
also a very recent preprint [84] for calculations of DOS 
in SrRuOa within an approach similar to ours. 

Figure 7 also presents for comparison the many-body 
density of states calculated for the material CaVOa 
[78], in which the correlations are believed to arise from 
Mott physics (Hund’s metal physics requires a higher oc¬ 
cupancy of the d level). The spectrum of the Mott ma¬ 
terial consists of a central peak corresponding to a uni¬ 
formly renormalized DFT DOS and clearly visible Hub¬ 
bard satellites at —1.5 and -1-3.5 eV. In contrast the spec¬ 
trum of the Hund’s material CaRuOa lacks clear Hub¬ 
bard sidebands, but as noted above does exhibit a more 
strongly renormalized and more highly structured quasi¬ 
particle peak. Although the Hubbard bands are more 
visible in the Mott material, the correlations are in an 
important sense weaker, with the mass enhancement of 
CaVOs only about 2-3 [78] in contrast to the ~ 7 that 
we find for CaRuOa. 


(a) Magnetization M and m* /m of SrRuO, 



temperature (K) 



temperature (K) 


FIG. 8. (Color online) (a) Solid curves: temperature depen¬ 
dence of orbitally-averaged spin-resolved mass enhancement 
of SrRuOs calculated at 1/ = 2.3 eV and J = 0.35 eV (left 
y axis). Dashed curve: magnetic moment M = (m^) calcu¬ 
lated for same parameters (right y axis). Inset: temperature 
dependence of squared magnetic moment = {{rtiz))^. (b) 
Orbitally-averaged mass enhancement vs. T for CaRuOa cal¬ 
culated at J = 0.35 eV and different U values; note CaRuOa is 
paramagnetic at all T. 


V. MAGNETIC PHASE OF SrRuOg 

In this section we show that the correlated Hund’s 
metal picture provides an adequate description of the 
magnetic phase of SrRuOa , preserving the successful de¬ 
scription of the magnetic properties obtained from DFT 
calculations [35, 43] while simultaneously providing a 
good account of the quasiparticle mass. While a restric¬ 
tion of the correlated subspace to the Ru t 2 g orbitals 
is adequate for most purposes, calculations within this 
scheme yield half-metals, with magnetic moments that 
saturate ai 2 ys per Ru site and with too-small values 
of the specific heat coefficient. We therefore treat a five- 
band model that includes both Cg and t 2 g-symmetry or¬ 
bitals, with the Cg states and eg-t 2 g interactions treated 
in a mean field approximation. [The intra-f 2 g interac¬ 
tions remain as specihed in Eq. (2)]. Inclusion of the Cg 
orbitals, which lie well above the Fermi level and play no 


important role in the calculation of dynamical quantities 
in the paramagnetic phase, acts to stabilize the magne¬ 
tization at partial polarization. 

The two panels of Fig. 8 compare the mass enhance¬ 
ments of SrRuOa and CaRuOa computed within the five- 
band model. The lower panel shows the T-dependence 
of the CaRuOa mass enhancements for three U values. 
The behavior is in agreement with the three-band com¬ 
putations discussed above: the mass increases as T is 
decreased, and below a D-dependent scale saturates at 
a C/-dependent value. At the J = 0.35 eV in this figure, 
we see that the mass is too small at U = 2.0 eV and ex¬ 
trapolates to a too-large value at U = 3eV, suggesting 
that a U « 2.3 eV provides a reasonable description of 
the physics. 

The upper panel shows the temperature dependence 
of the spin-resolved mass enhancement of SrRuOa for the 
parameters U = 2.3, J = 0.35 eV that provide a good de- 
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scription of CaRuOa , along with the temperature depen¬ 
dence of the calculated magnetic moment. The moment 
saturates to a value « 1.64Experiments report val¬ 
ues of the magnetic ranging from 0.8 to 1.6 /is per Ru site 
[77, 85, 86 ], with more recent experiments converging on 
values between 1.4 and 1.6 [87, 88 ]. DFT calculations 
[41, 42, 51, 69, 89] report magnetic moments consistent 
with experiment. We see that the Hund’s metal picture 
provides a similarly good level of agreement. The inset 
reveals that the magnetization has a mean-field-like tem¬ 
perature dependence linear in — T near the tran¬ 
sition; extrapolation to = 0 indicates a Curie tem¬ 
perature « 200 K slightly higher than the experimental 
Tc « 160 K. Fluctuation effects in a three dimensional 
magnet typically reduce the transition temperature by 
~ 30% relative to the DMFT value (see, e.g. Refs. 90 
and 91) so this value also is very reasonable. 

At temperatures above the magnetic phase transition 
the mass enhancement of SrRuOa is very close to that 
of CaRuOa. (The differences discussed above between 
the paramagnetic phase mass enhancements of the two 
materials become manifest only at low T where the mass 
enhancement is large and very sensitive to the spin dy¬ 
namics controlled by J and the density of states). As T is 
decreased below Tc we see that the mass enhancement in 
SrRuOs becomes spin dependent, taking different values 
for the majority and minority spin channels. The smaller 
value of the mass in the minority-spin channel may be un¬ 
derstood as a phase space effect. The dominant part of 
the interaction is between opposite spin species, embod¬ 
ied in the 5+5'“, part of the local interaction. Because 
we are dealing with a greater than half-filled band, the 
phase space available for a minority spin electron to scat¬ 
ter into a majority spin state is less than the phase space 
available for a majority-spin electron to scatter into a mi¬ 
nority spin state. More importantly, as the amplitude of 
the magnetic moment increases we see that the increase 
in mass is cut off, so the concavity of the m*/m{T) curve 
changes and the T dependence of the mass saturates be¬ 
low T ~ 50 K. This behavior is a natural consequence 
of the Hund’s metal physics, in which the large mass en¬ 
hancement arises from slow fluctuations of spontaneously 
generated local moments whose formation and dynamics 
is very sensitive to J and T [19, 23]. The quenching of 
these moments in the ordered phase then cuts off the in¬ 
crease of the mass enhancement. The values obtained for 
the majority-spin mass are in reasonable correspondence 
with experiment [69] although the contribution of the mi¬ 
nority spin channel to the overall specific heat requires 
further investigation. 

In Fig. 9, we present an estimate for the specific heat 
coefficient 7 = limT->.o C/T for both materials. We ob¬ 
tain the specific heat by using the quasiparticle approxi¬ 
mation Scr(w) ~ -I- (1 — Z~^)uj in the formula 


7 = 


3 


(Pk 

(^ 


Tr 


ImGcr(fc, zO )Z^^ 


( 4 ) 


which is obtained from the quasiparticle approxima¬ 



FIG. 9. (Color online) Specific heat coefficient 7 (Sommerfeld 
constant) calculated for SrRuOs as a function of temperature 
T at J — 0.35 eV and U values indicated. Inset: 7 for CaRuOa 
at U = 2.3eV and J = 0.35 eV. 


tion to the free energy F = —TrlnG. Here the spin- 
dependent quantities Sct,Go- and are 12 x 12 ma¬ 
trices, accounting for the twelve bands (three for each of 
the four Ru atoms in the unit cell.) Use of the quasi¬ 
particle approximation in Eq. (4) is strictly correct only 
in the low-T Fermi liquid regime, but is expected to be 
of the correct order of magnitude even at higher tem¬ 
peratures. In the orbital basis that minimizes the hy¬ 
bridization function (see Sec. H), Zg- is diagonal in site 
and orbital indices while the four Ru atoms are equiva¬ 
lent up to a rotation in orbital space. In this case, Eq. (4) 
simplifies to 


7 = 


ImG™™(r = 0,z0-) ^ 


3 ^ 

um 


( 5 ) 


where the sum is over spin and the d orbitals of one 
particular Ru atom. Figure 9 is produced using Eq. (5) 
with Zcrm and = ReEcrm(zw„ —>■ zO’*') obtained from 
DMFT. 

We see from Fig. 9 that the value U = 2.3 eV selected 
on the basis of the phase diagram and CaRuOa mass 
reproduces well the low-T experimental results (7 ~ 
30 mJ/molK^ for SrRuOa and 75 mJ/molK^ for 
CaRuOa [69, 77]). Larger (3.0eV) or smaller (2.0eV) 
value of U are inconsistent with experiment. Essential in 
obtaining this level of agreement is obtaining a correct 
estimate for the magnitude of the saturation moment. A 
too small moment {U = 2.0 eV) leaves active spin fluctu¬ 
ations which cause a further increase in the mass, while 
if the moment is too near to saturation {U = 3.0 eV) the 
7 is too small. 
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VI. CONCLUSIONS 

In this work, we have investigated the interplay of 
electronic correlations and lattice distortions in the per- 
ovskite ruthenates SrRuOa and CaRuOs using a den¬ 
sity functional treatment of the basic electronic states 
and treating correlations via dynamical mean-field theory 
with a CT-HYB impurity solver [63]. We determined the 
behavior of a general class of models motivated by the 
physics of the ruthenates, finding that ferromagnetism 
depends on (1) how far materials are from Mott insulat¬ 
ing phase and (2) positions of DOS peaks with respect 
to the Fermi level. The latter factor is more important 
for small and intermediate correlations, while the former 
controls the behavior at strong correlation. 

Our main results are presented in the phase diagram 
shown in Fig. 2. The choice of 17 = 2.3 eV value 
(far from critical value for the Mott insulating phase) 
and J = 0.35 eV gives a calculated mass enhancement 
for CaRuOs in reasonable agreement with experiment. 
The mass enhancement is a sensitive function of the 
distance of the material from the ferromagnetic phase 
boundary. Several experimental works propose that very 
weak disorder may induce spin-glass behavior in CaRuOa 
[77, 92, 93], suggesting that the material is very close to 
a magnetic phase boundary, consistent with our results. 
The same calculations predict that SrRuOa becomes fer- 
romagnetically ordered below a Curie temperature some¬ 
what less than room temperature. The onset of ferro¬ 
magnetism cuts off the low-T increase of the mass, and 
the calculated ferromagnetic-state value is found to be 
much smaller than that of paramagnetic CaRuOa and 
also in good agreement with experiment. We therefore 
concluded that the ruthenates are far from the Mott in¬ 
sulating phase and may be identified as Hund’s metals. 

A test of this picture is that the predicted density- 
of-states for CaRu03 is typical of a Hund’s metal: it 
shows only weak signatures of Hubbard bands despite 
the large renormalizations at low frequency; in contrast, 
the density of states of a material located near the Mott 
critical point would exhibit a much larger redistribution 
of weight away from the Fermi level and would display 
clear Hubbard sidebands. Closely related to this point 
is the local susceptibility shown in the inset to Fig. 4 
which is very strongly enhanced relative to the band the¬ 
ory value, demonstrating that the large renormalizations 
come from strongly enhanced local spin fluctuations. A 
further consequence of the identification of the materials 
as Hund’s metals is that the reason for the occurrence of 
ferromagnetism in SrRuOa and its absence in CaRuOa is 
the difference in density of states of the two materials. 

Our study has certain limitations. First, it employs the 
single-site dynamical mean field approximation. While 
this captures many important aspects of local energetics 
and material trends, it is not necessarily quantitatively 
accurate for d = 3 dimensional materials. In particular, 
the DMFT theory includes dynamical effects of momen¬ 
tum averaged spin fluctuations but does not include all 


of the effect of quantum critical spin fluctuations which 
are important near magnetic transitions [94]. Spin fluc¬ 
tuations are known to have particularly important con¬ 
sequences in quasi-two-dimensional materials and have 
been extensively discussed in the context of the two- 
dimensional ruthenate material Sr 2 Ru 04 ; their quanti¬ 
tative importance in d = 3 is less [94] but the actual con¬ 
tribution to physical properties, both in the ruthenates, 
and more generally in the context of DMFT, is an im¬ 
portant open question which warrants further research. 

Second, we have used the “frontier orbital” approxi¬ 
mation in which correlations are applied to the d-derived 
near-Fermi-surface bands. However, the good qualitative 
and even semiquantitative agreement between our cal¬ 
culations and experiment justifies this approximation a 
posteriori. Our results thus unambiguously indicate that 
the perovskite ruthenates are in the Hund’s metal class 
of materials, with strong correlation effects driven by the 
J rather than the U term of the interaction. This finding 
resolves the tension between the successful DFT account 
of the magnetism and the thermodynamic, transport and 
optical results indicating strong correlations. 

Third, we have neglected spin-orbit coupling which 
splits the sixfold degenerate t 2 g manifold into fourfold 
and twofold degenerate manifolds and tends to reduce 
bandwidths, and also can lead to spin anisotropy [32] 
and to interesting entanglement effects e.g. on super¬ 
conducting wave functions [95] . If the spin-orbit induced 
splitting (as renormalized by interactions) is large enough 
to significantly change the structure of the t 2 g mani¬ 
fold, for example by fully separating the twofold- and 
fourfold-degenerate submanifolds, then the basic picture 
presented here would not apply. On the other hand, if the 
spin-orbit coupling is not too large, then the tendency of 
J to favor high-spin states in which electrons are more 
or less equally distributed over orbitals will tend to sup¬ 
press spin-orbit effects. The crossover as spin-orbit cou¬ 
pling is increased at fixed J will have something of the 
character of a high-spin/low-spin transition and is thus 
expected to be abrupt. DFT calculations [96] for related 
material (Sr 2 Ru 04 ) yield a spin-orbit splitting of at most 
~ 200 meV and with no significant alteration of the basic 
properties of the bands. This, and the fact that calcula¬ 
tions of the kind presented here reproduce well the prop¬ 
erties of both pseudocubic and layered ruthenates (see 
Refs. 20, 82, and 84) suggests that in fact spin-orbit cou¬ 
pling does not have a strong effect on the basic electronic 
properties of the ruthenates. The issue however warrants 
more thorough investigation. 

We conclude by indicating few directions for the fu¬ 
ture work. Our calculations show an intricate interplay 
between lattice structure and correlation effects, medi¬ 
ated by the lattice-induced changes in the near-Fermi- 
surface density of states. Films show a different pattern 
of rotation and tilt than do bulk materials, and in films 
the rotation angles may be manipulated by strain. Cal¬ 
culations of the strain dependence of the mass enhance¬ 
ment and magnetic moment in thin films of SrRuOa and 


12 


SrRuO^ _ CaRuQ.^ 



FIG. 10. (Color online) DFT band structure (solid curves, 
blue online) and the corresponding MLWF fits (dots, green 
online) for SrRuOa (left column) and CaRuOs (right column) 
using their corresponding experimental orthorhombic struc¬ 
ture [52, 53]. (a) and (b): MLWF fits for t 2 g subspace, (c) 
and (d): MLWF fitting for t 2 g-eg subspace. The band struc¬ 
tures are plotted along the k path T^T^Y^X^U 
where F = 000, T = TTTrTr, Y = ttttO, X = ffTrO and U = Tfmr 
in the first Brillouin zone of the orthorhombic structure. The 
horizontal dashed line marks the Fermi level. 


CaRuOa will be very interesting to perforin and com¬ 
pare to experiment [46, 97] and to previous LSDA re¬ 
sults [44, 45]. Studying theoretically the ruthenates 
within a wide energy window to include oxygen p and eg 
bands, which are both close to the frontier t 2 g bands, is an 
interesting subject, too. Finally, the formalism and phys¬ 
ical picture provided here constitute a potentially use¬ 
ful starting point for investigations of impurity-induced 
magnetism in CaRuOa. 
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Appendix A: Wannier Fits 


Figure 10 presents the band dispersion obtained from 
our density functional calculations in several high sym¬ 
metry directions for SrRuOa and GaRuOa. The energy 
window has to be adjusted to capture all the states belong 
to the subspace. For the t 2 g-only subspace [Fig. 10(a) 
and 10(b)], an energy window from —3 to leV is used. 
For the t 2 g-eg subspace [Fig. 10(c) and 10(d)], the energy 
window ranges from —3 to 6eV. 

The frontier orbitals are seen to be well separated from 
other bands at higher energy, and as a result the MLWF 
fits are adequate over most of the relevant energy range. 
However, at the bottom of the t 2 g bands, some overlap 
with the oxygen p bands occurs, especially in the case 
of SrRuOa. The band overlap occurs near the F (zone 
center) point of the Brillouin zone, at which orbital char¬ 
acters are well defined and p-d hybridization is minimal. 
We have verified for the density of states projected to 
atomic orbitals of SrRuOa (not shown), the Ru d charac¬ 
ter is nonzero until ^ — 2.5eV, indicating that the Wan¬ 
nier fitting is thus reasonable even in this region of band 
overlap. 

As noted in the main text, if only the t 2 g manifold is 
included, the Wannier functions produced by the WAN- 
NIEr90 code [58] are aligned with the local octahedral 
axes and the DMFT hybridization function can be con¬ 
structed directly from the projection of the Kohn-Sham 
Hamiltonian onto the Wannier basis. If however all five 
d orbitals are included (as in our magnetic calculations) 
the orbitals produced by the the wannier90 code are not 
properly aligned to the local symmetry axes, and must 
be rotated, in order to minimize the off-diagonal terms 
in the DMFT hybridization functions. We find that the 
desired rotation is the one that diagonalizes the site-local 
terms in the projection of the Kohn-Sham Hamiltonian 
onto the Wannier basis. 

Our convention for the orbitals whose self-energy is 
shown in Fig. 5 is as follows. The Pnma structure has 
three lattice vectors conventionally denoted a, b and c. 
The orbitals are labeled in terms of pseudocubic x, y, z 
directions defined as those closest to 


d = x-\-y, 

(Al) 

b = 2i, 

(A2) 

c = x-y, 

(A3) 


and Fig. 5 presents the diagonal components of the self¬ 
energy for the Ru ion at position (0, 0, 0.5) in the system 
defined by the Pnma a, &, and c lattice vectors 
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FIG. 11. (Color online) Left column: the evolution of the 
Curie temperature for CaRuOa and SrRuOa with J increased 
from 0.25 eV to 0.5 eV and U = 3eV kept fixed; the fitting 
lines (solid for CaRuOs and dashed lines for SrRuOa) and 
all data points are included. Right column: expanded plots 
for smaller range of temperature near 0. The dotted line at 
J = 0.4 eV is the fitting line for CaRuOs if data of two lowest 
temperatures are neglected. The vertical dashed line marks 
zero temperature. 


Appendix B: Criteria for determining the electronic 
phases and phase bonndaries 

1. Ferromagnetic-paramagnetic phase bonndary 

To locate the ferromagnetic-paramagnetic phase 
boundary we follow our previous work [67] and compute 
the inverse magnetic susceptibility as x~^ = h/m with 
m the calculated magnetization and h an applied field 
chosen to be small enough that the m{h) curve is linear 
(typically 0.01 eV, but can be smaller at low tempera¬ 
tures). We perform the calculation at several tempera¬ 
tures, fit the result to a straight line X“^(T) = A{T — Tc) 
and determine the phase according to whether the T = 0 
extrapolation x“^(T = 0) = —AT^ is positive or nega¬ 
tive. 

For each point on the U-J phase diagram, we use 
at least four different temperatures, typically T = 
0.2,0.14,0.1 and 0.05 eV, to determine the Curie temper¬ 
ature (for some points close to the phase boundary, we 
go to lower temperatures). Within the magnetic phase 
the inverse susceptibility x~^ is typically linear in this 
range of temperature, thus is easily obtained. The 
phase boundary is then specified by linear interpolation 
between the points of lowest positive and largest negative 
Te- 

Figure 11 demonstrates the approach presenting 


X~^(T) for different J values at a fixed U = 3eV. We 
see clearly that SrRuOs has a greater tendency to ferro¬ 
magnetism than CaRuOa , with the difference being more 
pronounced at higher J. We note, however, within the 
paramagnetic phase or in the magnetic phase very close 
to the phase boundary, x~^ starts to bend away from the 
high-T linear extrapolation at low temperature, as can be 
seen from a close examination of the x~^ for CaRuOs for 
J = 0.33 and 0.4 eV. For these data we have pushed the 
CaRuOs calculation to the lower temperatures T = 0.025 
and 0.0125 eV. The bending away from the Curie-Weiss 
curve is a signature of the onset of Fermi-liquid coherence 
and leads to uncertainty in specifying the magnetic phase 
boundary. For example, at J = 0.4 eV, the linear extrap¬ 
olation gives Tc = 0.0006 eV if the points at T = 0.025 
and 0.0125 eV are included (solid line), but if these two 
temperatures are excluded (as in most of our calcula¬ 
tions for building the phase diagram), we would obtain 
Tc = 0.0052 eV (the dotted line in “zoom in” panel of 
Fig. 11). Therefore we expect the error bar for Tc of 
about 0.0045 eV. Our U-J phase boundary contains sim¬ 
ilar uncertainties, but determining the precise error bars 
At/, A J requires heavy calculations to go to ultra-low 
temperatures. The uncertainties arising from the onset 
of coherence do not affect the qualitative conclusions of 
this paper. 


2. Metal-insulator phase boundary 

We define whether the material is insulating or metallic 
according to whether or not the electron spectral function 
(many body density of states) A{uj) vanishes at the Fermi 
level w = 0 as T —)■ 0. We determine A{uj = 0) from the 
imaginary part of the measured Matsubara Green’s func¬ 
tion G{iujn) such that ImG(lw„) —)■ —ttA^uj = 0) when 
LOn —>■ 0. In practice, we observe ImG(ia;„) at several 
lowest Matsubara frequencies: if it bends towards zero, 
the state is insulator, whereas it goes away from zero, it 
is metallic. By fixing J and gradually increasing U, the 
critical value Uc is determined if the ImG(ia;„) bending 
changes at low frequencies. (See [98] for details). 

In single-site dynamical mean-field theory, the metal- 
insulator phase boundary has a complicated structure 
at low T, with a line of first order transitions emerg¬ 
ing T — 0 and there is second order transition at inter¬ 
action values Uc 2 ,Jc 2 and terminating at a critical end¬ 
point Uci, Jci, with Uci typically 0.8 — 0.9C/c2, there ex¬ 
ists a temperature Tmit where Uci{Tmit) = Uc 2 {Tmit)- 
Above this temperature there is only a crossover from 
metallic to insulating state [62]. We start all of our cal¬ 
culations from a metallic initial condition and the true 
metal-insulator transition is at a slightly lower tempera¬ 
ture than the Tmit we find, and the Uc we determine is 
closer to Uc 2 than Ud- 
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